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Abstract 

A scaling hypothesis for the n-particle spectral densities of the 0(3) non¬ 
linear sigma-model is described. It states that for large particle numbers 
the n-particle spectral densities are “self-similar” in being basically rescaled 
copies of a universal shape function. This can be viewed as a 2-dimensional, 
but non-perturbative analogue of the KNO scaling in QCD. Promoted to a 
working hypothesis, it allows one to compute the two point functions at all 
energy or length scales. In addition, the values of two non-perturbative con¬ 
stants (needed for a parameter-free matching of the perturbative and the non- 
perturbative regime) are determined exactly. 
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1. Introduction: An efficient way to describe the two-point function of some local 
operator O in a relativistic QFT is in terms of a Kallen-Lehmann spectral representation. 
The spectral density p(p) of O can be viewed as a measure for the number of degrees of 
freedom coupling to O at energy p. It decomposes into a sum of -n-particle contributions 


P(+) = > 

n 


(i) 


where, depending on the local operator under consideration, some of the n-particle contri¬ 
butions may vanish on the grounds of internal quantum numbers. In a theory with a single 
mass scale m one has p^(p) ~ <5(p — m) and p( n )(p), n > 2 has support only above nm, 
i.e. above the n-particle threshold. Once p(p) is known, the various (Minkowski space 
or Euclidean) two-point functions of O can be computed as convolution integrals with 
an appropriate kernel carrying only kinematical information. The dynamical problem 
consists in computing the n-particle spectral densities p^ n \p) of O. 

Here we shall be concerned with massive 1+1 dimensional QFTs and in particular with 
the 0(3) nonlinear sigma (NLS) model. The four most interesting local operators in this 
model are: The spin held, the Noether current, the energy momentum (EM) tensor and 
the topological charge (TO) density. Their spectral densities can be grouped into two 
families p\ n \p), n > 1, l = 0,1 according to their isospin. For n even/odd the p^ 
are the spectral densities of the EM-tensor/TC-density, respectively; similarly p, for n 
even/odd are the spectral densities of the Current/Spin, respectively. The following pieces 
of information are available for these spectral densities: (i) For small particle numbers 
n the functions p\ n \p) can be computed exactly by means of the form factor approach. 
In Q this has been done up to 6 particles, (ii) For all particle numbers n the p —► oo 
asymptotics of the n-particle spectral densities is known, and is given by 
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where the constants are computable from the integrals of the lower particle spectral 
densities [|TJ. The constants A^ are rapidly increasing with n. This implies that the 
p —> oo asymptotics of the full spectral densities (1) cannot be computed by naively 
summing up the asymptotic expressions (|2|), which in fact would be divergent, (iii) The 
large p asymptotics of the full spectral densities can however be computed in renormalized 
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perturbation theory (PT).Q One finds for the leading behavior 
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Subleading terms can also be computed, but not all of the overall constants are accessible 
to PT. In particular Ai := /L spin is an unkown non-perturbative constant. In the case 
of the TC density A tov is fixed by PT but its relation to the non-perturbatively defined 
spectral sum (1) is not. Equivalently the matrix elements of the TC density between the 
vacuum and some multi-particle state are defined only up to an unknown non-perturbative 
constant A 0 . 


Missing pieces of information about the spectral densities are: (iv) One would like to 
be able to compute the full spectral densities for all p > 0, not only their large p asymp¬ 
totics. This would allow one to compute the two-point functions at all energy/length 
scales. In terms of the spectral resolution (1) this amounts to knowing all the n-particle 
contributions, not only those with n < Hq for which the computation can be done explic¬ 
itly. (v) One would like to know the (exact) values of the non-perturbative constants Ao 
and Ai. Knowledge of these constants would allow one to match non-perturbative and 
perturbative information unambiguously. In this respect their role is similar to that of 
the mj A ratio 0. 


The purpose of this letter is to bridge the gap between the perturbative and the non- 
perturbative regime and to provide the missing pieces of information (iv) and (v). It is 
based on a remarkable self-similarity property of the n-particle spectral densities. For 
large n they appear to be basically rescaled copies of a “universal shape function” Yj(z). 
Explicitly 
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where Mj 71 - 1 and /f”'* are certain scaling parameters to be specified later. In the following 
we shall first give a precise formulation of the scaling law (|4|) and recall some of the 
evidence presented for it in flfl. Then we shall promote it to a working hypothesis and 
show that it has the following consequences: The UV behavior is consistent with PT; 
in particular those coefficients A° in (|3|) accessible to PT are reproduced. The non- 

1 The correctness of PT in this model has been challenged in ||. To simplify the exposition we shall 
assume the validity of PT for the UV asymptotics throughout this letter. 
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perturbative constants Ao and Ai are determined exactly and in the normalization 
given by 


Ao — 
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Finally candidate results for the two-point functions at all energy or length scales are 
obtained. 


2. Formulation of the Hypothesis: With hindsight to the asymptotics (j^) let us intro¬ 
duce 

R^ n \x) := me x p[ n \me x ), l = 0,1. (6) 

Here l — 0,1 as before correspond to the EM tensor & TC density and Spin & Current 
series, respectively. The graphs of these functions are roughly ‘bell-shaped’: Starting 
from zero at x — Inn they are strictly increasing, reach a single maximum at some 
x = > Inn and then decrease monotonically for all x > ff n ' 1 . The position (f"' 1 of 

the maximum and its value = R ( l n> (ff n ^) are two important characteristics of the 

function, and hence of the spectral density. Defining 
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both the value and the position of the maximum are normalized to unity. Initially Y^^z) 
is defined for (lnn)/£; < z < oo; in order to have a common domain of definition we set 

Y^ n \z) = 0 for 0 < z < (In n)/^. The proposed behavior of the spectral densities is as 
follows: 


Scaling Hypothesis: 

(a) (Self-similarity) The functions Y^ n \z), n > 2 converge pointwise to a bounded 
function Yi(z). The sequence of k-th moments converges to the k-th moments of 
Yi(z) for k + l = 0,1, i.e. 

]™o Y i H ( z ) = Y i( z ) ’ 

f* OO f* OO 

Jim y dzz k Yp\z) — J dzz k Y t (z). 

(b) (Asymptotic scaling) The parameters and scale asymptotically according 
to powers of n, i.e. 

^ n) ~ n 1+ai , M/ n) ~ Mi n~ 11 . 
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Feature (a) in particular means that for sufficiently large n the graphs of two subsequent 
members Yj n ^(z) and Y l < ' n \z) should become practically indistinguishable. This appears 
to be satisfied remarkably well even for small n = 4, 5, 6, as is illustrated in Figure [I] for 
the l = 1 series. 



Figure 1: Illustration of the self-similarity property of the rescaled l = 1 spectral densities. 
The plots show Y-j; n \z) (dashed) compared with (z) (solid) for n = 2,3,4, 5. 

The analysis of part (b) of the scaling hypothesis is more involved. It turns out that 
all but one of the exponents in (b) are fixed by self-consistency, and only this one has to 
be determined by fitting against the n < 6 particle data. The result is [[I 

7 i = 1 , «o = ai =: a , 

7o = 3 + 2a , a ~ 0.273 . (8) 

3. Consequences of the Hypothesis: Let no be the maximal particle number for which 
the spectral densities have been computed explicitly (at present n 0 = 6). Then (^) gives 
candidate expressions for all n > n 0 particle spectral densities so that one can evaluate 
their sum. Concerning the UV behavior of the sum notice that a finite number of terms 
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in the sum, each decaying according to (j2|), can never produce the different UV behavior 
in (|3|). However the infinite sum does. What is happening is that the partial sums 
(ln/r) 2_2i X^ 0 +i Pi\tA develop a plateau, i.e. are practically constant in a large interval 
hmm ~ /x ~ Hmax(N). When N is increased the plateau is prolonged and eventually 
reaches out to infinity, i.e. Hmax(N) —> oo for N —> oo (while is basically N- 

independent). The value of the plateau determines the asymptotic constants in (|3|). In 
those cases where the coefficients are accessible to PT the perturbative value is reproduced 
with an accuracy better than 1%. In addition one obtains the following two exact relations 
among the four constants 
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These relations reflect the linking between the even and the odd particle members of 
an isospin family, which results from the clustering relations obeyed by the exact form 
factors JT|. Since ( A CUII ) PT = 1/37r the first eq. in (|9|) gives H spm = Ai as in (||), while the 
second eq. fixes the physical normalization of the TC matrix elements in terms of that of 
the EM tensor, which in our conventions amounts to Aq = 1/4. 



m\x 


Figure 2: Two-point function of the topological charge density. Comparison MC data 
and form factor result. Test of Ao = 1/4. 

The Ao value can be tested independently by means of MC simulations. In Figure |2] the 
results for the two-point functions of the TC density is shown - once computed via the 
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form factor approach, with the absolute normalization fixed according to (|5j) and once via 
MC simulations. The simulations were done using the cluster algorithm of [}3j with the 
standard action and the geometrical definition of the TC density. The data in Figure 2 
correspond to a 460 square lattice at inverse coupling f3 = 1.80 (correlation length 65.05). 
The statistical errors are smaller than the size of the dots. The nice agreement confirms 
A 0 = 1/4 and hence provides further support for the scaling hypothesis. 

The plateau-phenomenon for the spectral densities has a counterpart for the two-point 
functions. In momentum space the latter can be computed as Sticltjes transforms of 
the corresponding spectral densities. Inserting the decomposition m one obtains for the 
n < N particle approximants 


V r co 
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For n > no the n-particle spectral densities are evaluated by means of the scaling hypoth¬ 
esis. The N —> oo limit of (10) yields a candidate for the exact two-point function. 



y = In p/m 


Figure 3: Spin two-point function: Ratio of In(u) approximants and 2-loop perturbation 
theory Ipr(y ) versus y = In p/m. 
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For the case of the spin held the result is shown in Figure in units of the 2-loop PT 
result for the same quantity. The ‘exact’ N = oo curve approaches its asymptotic value 
(here normalized to unity by using Ai = 4/(37r 2 ) for the normalization of the PT result) 
fairly slowly. Demanding a 1% accuracy for PT one finds that there are two such regimes, 
one at an intermediate energy scale 4 ~ y ~ 9 and a second one for y ~ 75. A similar 
behavior is found for the current two-point function. Thus, somewhat surprisingly, PT 
does not necessarily improve monotonically with increasing energy. As a reservation one 
must add however that these are 1% — 2% effects which may be affected by subleading 
terms in the scaling law (b). (Inclusion of a 3-loop PT contribution does not remove the 
effect; sublcading terms arising from extrapolating Y^ n+1 \z) _ Y^ n \z) have been taken 
into account in Figure 3.) 

4- Relation to KNO Scaling: A bonus of the above scaling hypothesis is that it implies 
KNO-like scaling laws for multi-particle production processes (but not vice versa). Of 
course there is no particle production in an integrable model, but we can enlarge our 
model-world by a “weak” sector so that a general state will be of the form |s;w), where 
s is the “strong” (sigma-model) part of the state and w is the “weak” part. Adding a 
current-current interaction term, production processes \s m ,w) —► \s']w') (where the states 
s and s' have different particle numbers) become possible. The corresponding transition 
amplitude, to lowest order in the new interaction, reads 

a = j d 2 xr(x)( s ' b»| s ), (ii) 

where the Fourier transform of Z“(x) has support at the momentum Q , the “weak” mo¬ 
mentum transfer. 

The simplest production process is the two-dimensional analogue of the e + e _ annihila¬ 
tion. We model this process by choosing Q 2 > 0 and s = |0). Summing over all (discrete 
and continuous) quantum numbers of the n-particle state s', the probability distribution 
for producing n particles at energy y = ^[Q 1 becomes independent of the “weak” part of 
the process and is proportional to the current spectral density. The proportionality factor 
involves the coupling constant of the external current-current interaction, which drops out 
when considering the conditional probability P 2n for having exactly 2 n particles produced, 
once the process has taken place at all. One has 

OO 

P 2 n=-df n) (/i), E P 2n = l, (12) 

K n=l 
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where the second eq. fixes the normalization constant k. Similarly production processes 
upon perturbation with the EM-tensor can be studied, in which case (/i) appears in 


flT2|). Using our scaling hypothesis the energy dependence of the probability distribution 
can in both cases be expressed, for large p (where one can approximate the sums by 
integrals), in terms of n = 2 uP 2 n ~ (ln//) 1 /( 1+a \ the average number of particles 

produced. The asymptotic distribution takes the KNO-scaling form || 


nP 2 n = f 


n 


n 


(13) 


The scaling function f(q) is given in terms of the universal shape function as 




hi 


1+aN 


,2 hi q, 


(14) 


where the parameters hi and hi are the (yy^ — l) th and (yyy — l) th moments of the 
universal shape function l)(z), respectively, and the exponents are given in (H). The case 
l = 1 corresponds to the current perturbation and l = 0 to the perturbation by the EM- 
tensor. These KNO-type scaling laws, however, are only valid for simultaneously large 
particle numbers and large energies. In contrast, the scaling hypothesis (a), (b) for the 
spectral densities is valid for all energies, and in particular is non-perturbative in nature. 
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